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The electronic structure of a graphene superlattice composed by two periodic regions with different Fermi velocity, energy gap 
and electrostatic potential is investigated by using an effective Dirac-like Hamiltonian. It must be expected that the change of 
the Fermi velocity in one region of the graphene superlattice is equivalent to changing the width of this region keeping the Fermi 
velocity unchanged, provided that the time taken to charge carriers cross the region is the same. However, it is shown here that these 
two systems are not equivalent. We found extra Dirac points induced by the periodic potential and their location in the k space. It 
is shown that the Fermi velocity modulation breaks the symmetry between the electron and hole minibands and that it is possible to 
control the behavior of the extra Dirac points. The results obtained here can be used in the fabrication of graphene-based electronic 
devices. 
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1. Introduction 

Since its first successful experimental realization 03], 
graphene, an one-atom thick layer of carbon atoms arranged 
in a honeycomb lattice, has attracted a great deal of attention. 
Such interest is due to its unusual physical properties and po¬ 
tential application J1S1. The charge carriers in graphene, 
for instance, have an extremely high mobility and can be easily 
controlled by applying a gate voltage, which makes graphene 
a promising material for the fabrication of electronic devices. 
Controlling the energy band structure of graphene can help its 
application in electronics. It can be done, for instance, by ap¬ 
plying and external periodic potential, i. e., a superlattice. 

The periodic potential structure in graphene can be gen¬ 
erated by different methods, such as electrostatic potentials 
HSHiSHmdEi, magnetic barriers HEEH 
or the combination of both II18L 11911 . Despite the difficulty of 
fabricating graphene under nanoscale periodic potentials, it was 
already realized experimentally 0 eh mm 00 . Al¬ 
though significant advances have been made in the understand¬ 
ing of how a periodic potential influences the electronic prop¬ 
erties of graphene, very little further attention appears to have 
been given to the effects of a periodic modulation of the Fermi 
velocity in graphene supeiiattices l2fl[ 27j 28. 29. 301. 


In our previous work 0, we have found that it is possible 
to control the energy gap of graphene with Fermi velocity engi¬ 
neering when there is an energy gap modulation. It happens be¬ 
cause the change of the Fermi velocity in a region of graphene 
is equivalent to changing the width of this region and keeping 
the Fermi velocity unchanged, provided that the time taken to 
charge carriers cross the region is the same. So, at the same way 
that it is possible to tune the energy gap of graphene controlling 
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the width of two regions with different energy gap, it is possi¬ 
ble to do the same changing the Fermi velocity. However, as 
will be shown here, this equivalence is related only with the en¬ 
ergy gap, and does not remain for other electronic properties of 
graphene. In this paper we investigate the electronic structure of 
a graphene superlattice with a periodic modulation of the Fermi 
velocity and energy gap. Including the periodic potential, extra 
Dirac points appear. We show that the Fermi velocity affects the 
extra Dirac points in a different way of changing the widths of 
the regions, being possible to control the behavior of the extra 
Dirac points in different ways. 

The paper is organized as follows. In Sec. 2 we solve 
the Dirac-like equation for graphne with a position dependent 
Fermi velocity, energy gap and electrostatic potential and find 
the dispersion relation. We also discuss how to obtain such a 
system experimentally. In Sec. 3 we analyze the electronic 
structure in two different cases. First we consider the Fermi 
velocity constant and investigate the electronic structure with 
a unequal well and barrier widths. Then, we consider a Fermi 
velocity modulation and an equal well and barrier widths. We 
compare these two cases and show that they are not equivalent. 
In Sec. 4 the paper is summarized and concluded. 

2. Model 

The effective two-dimentional Dirac Hamiltonian for a 
graphene superlattice with a position dependent energy gap and 
Fermi velocity is written as 

H = -ih(s]v F (x)o- x d x 4 v f (x) + v F (x)cr y d y ') 

+V(x)l + A(x)cr z , (1) 

where V(x) is an external periodic potential, A(x) is the energy 
gap, <x, are the Pauli matrices acting on the pseudospin related 
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to the two graphene sublattices, 1 is the 2x2 unitary matrix 
and v(x ) is the Fermi velocity. The first term on the Hamilto¬ 
nian above is modified in relation to the usual Dirac operator 
for graphene in order to have a Hermitian operator 13111 . We 
are considering that the Fermi velocity, the energy gap and the 
potential change only in the x direction. A schematic diagram 
is shown in Fig. Q] The graphene is deposited on a heterostruc- 
tured substrate composed by two different materials, which can 
open different energy gaps in different regions of the graphene 
sheet that will be denoted by Ai =0 and A 2 = A. The modula¬ 
tion of the Fermi velocity can be obtained in graphene by plac¬ 
ing metallic planes close to the graphene sheet, which will turn 
electron-electron interactions weaker and, consequently, mod¬ 
ify the Fermi velocity [32, 3- The Fermi velocity in each 
region will be denoted by vi and \ ! 2 - The electrostatic potential 
for the two regions will be denoted by Vi — 0 and V3 = V. The 
period of the system is a + b. 



Figure 1: Schematic diagram of a graphene sheet deposited on a heterostruc- 
tured substrate composed of two different materials which can open different 
energy gaps in different regions of the graphene sheet inducing a periodic mod¬ 
ulation of the energy gap in graphene. Metallic planes are placed close to 
graphene sheet which induce a periodic velocity barrier. The period of the 
graphene superlattice is a + b. 


Following the calculations done in 0, is straightforward to 
obtain that the dispersion relation is given by 


cos (k x l) 


■ 


cos(k\a) cos(k 2 b) 


k 2 h 2 viv 2 + E(V - E) 
H 2 V\V2k\k2 


sin(kia) sinfefr) , (2) 


where k! = (E 2 /h 2 v 2 -k 2 ) 112 , k 2 = ([(V-E) 2 -A 2 ]/h 2 v 2 2 ~k 2 ) l/2 
and we have defined l = a + b. Note that at V — A = 0 and 
Vi = v'o = \’h we recover the linear dispersion relation of a 
graphene sheet. 


3. Electronic Structure 

For the sake of comparison, in this section we will investi¬ 
gate the electronic structure of two special cases. First, we will 
consider a constant Fermi velocity and analyze the effects of 
an unequal well and barrier widths. Then, we will consider an 
equal well and barrier widths and see the influence of the Fermi 
velocity modulation on the electronic structure of the graphene 
superlattice. In what follows, we will consider a constant period 


of the graphene superlattice a+b - I = 60 nm and we shall con¬ 
centrate our discussion on the valence and conductance mini¬ 
bands only, assuming the Fermi level to be in between at any 
value of V. 

3.1. Unequal barrier and well widths 

In this section we will analyze the effects of an unequal well 
and barrier widths and we will consider v\ = V 2 - In Fig. [2] 
(a) the electron and hole minibands are plotted at a — b — 30 
nm (black), a — 20 nm and b - 40 nm (blue) and a = 40 nm 
and b = 20 nm (red) with A = 15 meV, k y = 0 and V - 0 
(continuum lines) and V 4 0 (dashed lines). One can see that, 
at V = 0, it is possible to tune the energy gap of graphene by 
changing the well and barrier widths. Increasing the potential 
the electron and hole minibands are shifted up, but not equally, 
which changes the energy minigap that may be zero. It is clear 
in Fig. [2 ](i>), where the electron and hole minibands are plotted 
as a function of V with k x = k y = 0. It is possible to see that 
the energy minigap oscillates and it is equal to zero at discrete 
values of V, which depends on the well and barrier widths. In 
the dashed lines in Fig. [2] (<a) we consider the first value of V 
that closes the minigap for the three cases. 



Figure 2: The dispersion relation 0 with ci — b — 30 nm (black), a = 20 nm 
and b = 40 nm (blue) and a = 40 nm and b = 20 nm (red) with vi = V 2 = 10 6 
m/s, A = 15 meV and k y = 0. ( a ) The electron and hole energies as a function 
of k x , with V = 0 (continuum lines) and V = 139.4 meV (dashed black lines), 
V = 157.14 meV (dashed blue lines) and V = 156 meV (dashed red lines), 
revealing the possibility of closing the energy gap with the potential. ( b ) The 
electron and hole minibands as a function of V with k x = k y = 0, which shows 
the oscillation of the energy gap. 


In Fig. [3] the electron and hole minibands are plotted as a 
function of k y for k x = 0. The colors represent the same cases 
that in Fig. [2] In order to be clear the influence of the periodic 
potential and the unequal well and barrier widths, the energies 
in this figure are counted from the contact point or from the 
minigap center for all values of V. In Fig. 0(a) we consider the 
first value of V that closes the minigap for the three cases. One 
can see that the contact point in k y direction has a parabolic be¬ 
havior, while in the k x direction it is conical. So, the minibands 
have a lens-like shape. It can be seen also that when the well 
and barrier widths are equal, the electron and hole minibands 
are symmetric related to £ = 0. However, this symmetry is 
broken when the well and barrier widths are different. 

In Fig. 0(b) we consider V = 207.89 meV (black), V = 
234.12 meV (blue) and V = 233.27 meV (red), which are val¬ 
ues between the first and second value of V that close the gap at 
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ky = 0. One can see that even though an energy minigap opens 
at k y = 0, extra Dirac points appear at k y 4 0. These extra 
Dirac points appear when V exceeds a critical value V c , which 
is the first value that closes the minigap, and never disappear. 
So, from V = V c , the graphene becomes gapless. One can note 
that the extra Dirac points that appears at k y 4 0 are not at the 
Fermi level when a 4 b. When a > b (a < b) the Dirac points 
are shifted up (down) the Fermi level. 



Figure 3: The dispersion relation 0 as a function of k y , with k x = 0. The colors 
here represent the same cases that in Fig. [2] We consider (a) the same values 
of V that in the dashed lines in Fig. Qand (b) V = 207.89 meV (blabk lines), 
V = 234.12 meV (blue lines) and V = 233.27 meV (red lines). 


In order to localize the contact points in k space let us take 
advantage of the implicit function theorem. One can conclude 
that at the contact points, where there is an intersection of the 
bands, the gradient (Jacobian) of the dispersion relation <[2j> 
should be zero. It happens only if sinAqa = sin k 2 b = 0 and 
cos k\ci = cos k 2 b = +1. So, one can write 


k\a = 


' E 2 

, fr 2 v 2 F 



a = nn 


(3) 


and 


k 2 b 


'(V - E) 2 - A 2 



b = nn, 


(4) 


where n is an integer different from zero. Subtracting ([4} from 
0, one gets 


_ V 2 - A 2 n 2 n 2 h 2 v 2 F / 1 1 \ 

~ 2V + 2V \a 2 ~b 2 )’ 


(5) 


which gives the values of E where there are contact points. Re¬ 
placing the equation above in Eq. 0 one obtains 


In order to be clear the location of the contact points it is 
important to explain that when V exceeds V\, two extra Dirac 
points appear at k y 4 0. When V = V 2 the energy minigap 
at k y = 0 closes and the system has now three contact points. 
When V exceeds V 2 the minigap at k y = 0 opens but two extra 
Dirac points appear, leaving the system with four contact points 
and so on. When V = V„ or V„ < V < V n+ \, E n and k yn give the 
location of the contact point nearest to k y = 0 whereas E\ and 
k yi give the location of the contact point farthest to k y = 0. So, 
the location of the contact points are 

(£, ky, ky) = (E n , 0, kyj, (£„_ 1,0, kyyj, ..., (E U 0, kyy) . (8) 
3.2. Periodic velocity barrier 

In this section we will consider a periodic velocity barrier 
with a = b = 30 nm. In Fig. [4] we consider iq = v 2 = 10 6 
m/s (black), Vi = 1.5 X 10 6 m/s and v 2 = 0.75 x 10 6 m/s (blue) 
and tq = 0.75 x 10 6 m/s and v 2 = 1.5 x 10 6 m/s (red) with 
A = 15 meV and k y = 0. We recover the result obtained in 
0, which showed that it is possible to tune the energy gap of 
graphene by Fermi velocity engineering. One can note that the 
Kg. a is exactly the same as Fig. [2] It happens because the 
values for the Fermi velocity were chosen here in such a way 
that the time taken to charge carriers cross the regions of the 
graphene superlattice is the same that in the case shown in Fig. 
[2] So, if we consider an one-dimensional case, the periodic 
velocity barrier is equivalent to the unequal well and barrier 
widths problem. 



Figure 4: The dispersion relation (2) with vi = V2 = 10 6 m/s (black), vi = 
1.5 x 10 6 m/s and V2 = 0.75 x 10 6 m/s (blue) and vi = 0.75 x 10 6 m/s and 
V2 = 1.5 x 10 6 m/s (red) with a = b = 30 nm, A = 15 meV and k y = 0. ( a ) The 
electron and hole energies as a function of k x , with V = 0 (continuum lines) and 
the dashed lines have the same values of V that in Fig. [2] ( b) The electron and 
hole minibands as a function of V with k x = k y = 0, which shows the oscillation 
of the energy gap. 


ky f 


El- A 2 

h 2 v 2 


( 6 ) 


As all contact points are located at k x = 0, with the equa¬ 
tions above we have the complete location of the contact points. 
From the zeros of Eq. 0 one obtains, 

T7 nnhvp (nnhv\ 2 

K =— + i(—) +A - <7) 

which are the values of V that close the minigap at k y = 0. The 
critical value of V is V c = V\. 


However, with k y 4 0 the two cases discussed here are not 
equivalent. In Fig. [5]we plotted the electron and hole minibands 
as a function of k y with k x - 0, where the colors represent the 
same cases that in Fig. |4] In Fig. 0(a) we consider the first 
value of V that closes the gap. One can see that when vq 4 v 2 
the electron and hole minibands are asymmetric, one becoming 
narrower than the other. When one exceeds this value of V, an 
energy minigap opens at k y = 0 and extra Dirac points appear at 
k y 4 0, as shown in Fig. 0(/?). But, in contrast to the extra Dirac 
points that appears in the case analyzed in the last section, the 
extra Dirac points here are all located at the Fermi level. Thus, 
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the periodic Fermi velocity barrier can not shift up or down the 
extra Dirac points, as can be done by changing the well and 
barrier widths. 



Figure 5: The dispersion relation 0 as a function of k y , with k x = 0. The colors 
here represent the same cases that in Fig. \4\ We consider (a) the same values 
of V that in the dashed lines in Fig. [2] and ( b ) the same values of V that in Fig. 
0 ( 6 ). 


The location of the extra contact points can be found at the 
same way as in the last section. One can note that k\ = k 2 when 


Vv 2 - vi Jv 2 v 2 - A 2 v 2 + A 2 Vj 

E = E 0 =- y ~ , ---. (9) 

v- - v 2 

As all contact points are located at k x — 0, one can see that Eq. 
(0 with a = b,k x = 0 and E = Eq is given by 


1 = COS 2 (fcifl) + 


k 2 tv 1 V\V 2 + E 0 (v - E 0 ) 
h 2 v\v 2 k\ 


■ sin“(kia) , (10) 


which is satisfied when k\ = mn/a, where m is an integer dif¬ 
ferent of zero, because m = 0 implies k\ = 0, which makes 
the denominator h 2 v\v 2 k\ in Eq. ( TTTTb vanishes. This condition 
leads to 


ky — ky m — 



which gives the values of k y where the contact points are lo¬ 
cated. As it was already observed, in contrast to the case dis¬ 
cussed in the previous section, in this case for a given value 
of V all contact points are located at the same value of energy 
Eq. The zeros of the equation above give the contact points at 
k y = 0. So, one can write 


mnhvi + .Jm 2 n 2 h 2 v 2 + a 2 A 2 

V m = -----, (12) 

a 

where the critical value of V wherewith the graphene superlat¬ 
tice becomes gapless is V c - V\. 


4 . Conclusions 

We analyzed the electronic structure of a graphene superlat¬ 
tice with a piecewise constant periodic potential, energy gap 
and Fermi velocity. We consider that the periodic energy gap is 
generated by an appropriate substrate composed by two differ¬ 
ent materials that open different energy gaps in different regions 


of the graphene and the Fermi velocity is modulated by metallic 
planes placed close to the graphene sheet. 

We compared two cases, one with an unequal well and barrier 
widths keeping the Fermi velocity constant and other with a pe¬ 
riodic Fermi velocity barrier keeping the well and barrier widths 
equal. It was shown that with k y = 0 the two cases are equiva¬ 
lent, being possible to tune the energy gap of the graphene su¬ 
perlattice by changing the well and barrier widths or by Fermi 
velocity modulation. It was also found that the energy gap os¬ 
cillates when the potential increases continuously. However, 
with k y + 0 the two cases considered here are not equivalent. 
We found that extra Dirac points appears at k y + 0 when V ex¬ 
ceeds the first value Vi that closes the energy gap at k y = 0. 
As a consequence, the graphene superlattice becomes gapless 
after V reaches a critical value V, = V\. We showed that in 
both cases the electron and hole minibands become asymmet¬ 
ric. But the extra Dirac points for the case of unequal well and 
barrier widths are shifted up or down the Fermi level, which 
does not happen with a periodic Fermi velocity barrier. The 
results showed here provide different ways of controlling the 
behavior of the extra Dirac points in graphene and can be used 
in the fabrication of graphene-based devices. 
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